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Demonstrating the quantum-confined Stark effect (QCSE) in silicon nanocrystals (NCs) embedded 
in oxide lias been rather elusive, unlike the other materials. Here, the recent experimental data from 
ion-implanted Si NCs is unambiguously explained within the context of QCSE using an atomistic 
pseudopotential theory. This further reveals that the majority of the Stark shift comes from the 
valence states which undergo a level crossing that leads to a nonmonotonic radiative recombination 
behavior with respect to the applied field. The polarizability of embedded Si NCs including the 
excitonic effects is extracted over a diameter range of 2.5-6.5 nm, which displays a cubic scaling, 
a = cDnCi "^i^h c = 2.436 x lO"^'^ C/(Vm), where -Dnc is the NC diameter. Finally, based 
on intraband electroabsorption analysis, it is predicted that p-doped Si NCs will show substantial 
voltage tunability, whereas n-doped samples should be almost insensitive. Given the fact that bulk 
silicon lacks the linear electro-optic effect as being a centrosymmetric crystal, this may offer a viable 
alternative for electrical modulation using p-doped Si NCs. 

PACS numbers: 73.22.-f, 78.67.Bf, 71.70.Ej 



; I. INTRODUCTION 

The Stark effect has evolved within the previous cen- 
tury into a poviferful spectroscopic tool for the solidsii 
In the case of bulk semiconductors, shallow free excitons 
become easily ionized which limits the strength of the 
applied fields. A more robust variant of the Stark effect 
is the so-called quantum-confined Stark effect (QCSE) 

I where the carriers are trapped in a quantum well or a 
lower dimensional structure.^ In this regard, the quantum 
dots or nanocrystals (NCs) are preferred so as to take 
advantage of the full three-dimensional confinementi^ 
As a matter of fact, the electroabsorption studies were 
initiated quite early with CdS2:Sei_j; NCs embedded 
in a glass matrixi^ The QCSE activity in group-ITVI 
NCs^"— was soon extended to group-IITAs NCs/^i^ A re- 
lated noteworthy achievement was registering photolu- 
minescence (PL) from a single quantum dot within an 
ensemblefi^ followed by probing the QCSE from a single 

; dotJi 

On the technologically important front of group-IV 
materials, a recent breakthrough was the announcement 
of QCSE in germanium multiple quantum wells sand- 
wiched between SiGe barrier layers.''-^ The drawback of 
this structure is the small band offset of the barrier re- 
gions which limits the applied reverse bias before car- 
' rier tunneling sets in. Furthermore, it suffers from the 
polarization-dependent response discriminating between 
TE and TM polarizations; both of these shortcomings 
are inherently carried over to Si/Ge self-assembled quan- 
tum dotsii^ Si NCs embedded in oxide, not only offer 
remedy to both of these problems, but also due to its in- 
sulating host matrix it can withstand very high electric 
fields. Surprisingly, even though nanosilicon has become 
an established field/* the QCSE activity in this system 
has been quite overlooked. In some of the early electro- 
luminescence and photoluminescence studies on Si NCs, 
the precursors of QCSE was reported as a small redshift 



which was however taken over by a strong blueshift J^ii^ 
As another indirect measurement of QCSE in Si NCs, 
Lin et al. announced all nm redshift under strong illu- 
mination, but without an external bias, which they at- 
tributed to a build up of an internal electric field due 
to capture of carriers in NCs."'^^ Only very recently, the 
direct measurement of QCSE under an external field 
in Si NCs was achieved on ion-implanted samples that 
yielded as large as a 40 nm redshift at cryogenic tempera- 
tures, which remained to be easily detectable at the room 
temperature^^ This much delayed progress may never- 
theless become crucial for the electronically controllable 
silicon-based photonics and especially for optical modu- 
lators; the latter has been a real challenge, as bulk silicon, 
being a centrosymmetric crystal, lacks the Pockels effect 
which leaves the plasma effect as the main route for elec- 
trical modulation] Very recently, a GeSi electroab- 
sorption modulator has been announced that makes use 
of the bulk Franz-Keldysh effect of germanium enhanced 
under tensile strain. Amidst these developments, the 
present understanding on the QCSE and electroabsorp- 
tion in Si NCs remains to be quite insufficient so as to 
address whether it can offer a viable alternative to the 
existing and emerging ones. 

In this work, we aim for an assessment of these 
prospects from a rigorous atomistic point of view, start- 
ing with the recent QCSE experiment and extending our 
analysis to both fundamental as well as applied direc- 
tions. First, we theoretically show that, the high ly p ro- 
nounced luminescence shifts as measured in Ref. [l8| un- 
ambiguously originates from the Stark effect. In so do- 
ing, the importance of the excitonic effects is emphasized 
for larger NCs. The detailed explanation of the emis- 
sion strength as a function of Stark field reveals the in- 
tricate interplay of the single-particle Stark shifts, their 
level crossings and the electric field dependence of the 
oscillator strengths. From a fundamental point of view, 
in this context the most important physical quantity is 
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their polarizability. However, this is a subject which has 
not been discussed so far in the Uterature. Therefore, 
we provide the polarizabihty of embedded Si NCs for the 
useful 2.5-6.5 nm diameter range and furthermore, our 
exhaustive computations are expressed in simple expres- 
sions to enable their use by other researchers. Finally, we 
complete our theoretical analysis with the intraband elec- 
troabsorption properties of Si NCs under a strong electric 
field, where we predict the marked discrepancy between 
the n- and p-doped NCs. 

Our computational framework is a semiempirical 
pseudopotential-based atomistic Hamiltonian-- in con- 
junction with the linear combination of bulk bands 
(LCBB) as the expansion basisi^i^ The strong Stark 
field is included directly to the Hamiltonian without any 
perturbative approximation. This is the current state- 
of-the-art theory for this system size, which is far ad- 
vanced compared to effective mass and envelope func- 
tions approaches'^"— (for a critical account of the latter, 
see Ref. [28h and moreover is not amenable by ab ini- 
tio techniques)^ The competence of this technique has 
been well tested; in the case of embedded Si and Ge 
NCs, this has been employed to study interband and 
intraband optical absorptions^ the Auger recombina- 
tion, and carrier multiplication;^ furthermore, predic- 
tions for the third-order nonlinear optical properties us- 
ing the same theoretical approach^'^ have been indepen- 
dently verified experimentallyi^ The paper is organized 
as follows. A description of the theoretical method is 
given in Sec. II. The QCSE, polarizability, and intraband 
electroabsorption are analyzed in Sec. III. Main conclu- 
sions are presented in Sec. IV. 



II. THEORY 

In the LCBB approach^ the NC wave function with 
a state index j is expanded in terms of the bulk Bloch 
band (n) and the wave vector (k) as. 



the pseudopotential framework, given by 



n.k,fi 



(1) 



where N is the number of primitive cells within the 
computational supercell, C'^ - is the expansion coefH- 

n,k.J 

cient set to be determined, and fj, is the constituent bulk 
material label pointing to the NC core and embedding 
medium, -(r) is the cell-periodic part of the Bloch 

states which can be expanded in terms of the reciprocal- 
lattice vectors, {G} as 
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nk 



G 



JG-r 



(2) 



where is the volume of the primitive cell. The Hamil- 
tonian has the usual kinetic and the ionic potential parts, 
the latter for describing the atomistic environment within 



H ^T + Vpp 



2mn 



E Wi:{R,)v^^[r~R,-d^^) ,(3) 



where toq is the free electron mass, WI^{Rj) is the atomic 
identity coefficient that takes values or 1 depending on 
the type of atom at the position Rj — d^, here Rj is the 
primitive cell coordinate and is the displacement of 
this particular atom within the primitive cell, vf^ is the 
screened spherical pseudopotential of atom a of the ma- 
terial /i, the latter distinguishes the NC and the matrix 
regions. 

The formulation can be cast into the following gener- 
alized eigenvalue equation^^^^ 

EH ^ - r"^ - = E 9 - - r*^ ^ (A) 
n' k' fi' ,nkfi n^kj ^ / ^ n' k' fi' ,nkfi n,kj ' ^ 



n.k.fj, 

where. 
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the overlap part in the generalized eigenvalue equation is 
of the form 



S ,t, , r = in'k' ix'\nkix) 

n' k' ^' .nk^ \ ^ i ^ / 



The Si NC is intended to be embedded in silica, repre- 
sented by an artificial wide band gap host matrix that has 
the same band-edge line up and the dielectric constant, 
but otherwise lattice-matched with the diamond struc- 
ture of Si4i We refer to Ref. [Sll for the other technical 
details of the implementation of the electronic structure, 
including the form of the pseudopotentials for the NC 
and matrix media. 

For the study of QCSE, we treat the strong external 
field in the same level as the other terms of the atom- 
istic Hamiltonian (i.e., nonperturbatively). At variance 
with the electrostatic model used in Ref. [iSl, we assume 
that an individual spherical Si NC under consideration is 
embedded in a uniform medium having a constant per- 
mittivity for silicon rich oxide. This is justified by the 
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FIG. 1. (Color online) The electrostatic setting of the embed- 
ded NC under a dc external field. 



spatial distribution of the light-emitting centers in Si- 
implanted Si02'^ In connection to the actual samples, 
we inherently assume that the NCs are well separated 
which applies safely to NC volume filling factors of about 
10% or less. The basic electrostatic construction of the 
problem is presented with the assumption that the NCs 
are well separated in Fig. [TJ If we denote the uniform 
applied electric field in the matrix region asymptotically 
away from the NC as Fq, then the solution for electro- 
static potential is given in spherical coordinates by^ 

{— j^i^orcos6', r < a 

-Forcos0+ (f^)Fof^cos0, r>a ' 

. . . 

where e = eNc/cmatrix is the ratio of the permittivities 
of the inside and outside of the NCs. Hence, this ex- 
pression accounts for the surface polarization effects due 
to dielectric inhomogeneity which partially screens the 
external (i.e., dc Stark) field. The effect of this exter- 
nal field can be incorporated by adding the VJjxt = e$ 
term to the potential-energy matrix elements. A com- 
putationally convenient recipe in the context of QCSE 
is to assume that external potential is relatively smooth 
so that its Fourier transform can be taken to be band 
limited to the first Brillouin zone of the underlying unit 
cellj^ This results in the following LCBB matrix element 

(n'k'^^'\%^,\n^ = e^k - k') ^ B';^,. {g')* B';^-, (g) 

G,G' 

^^^^^bM(G + k-G' -k'), (6) 

here $(fc) is the Fourier transform of ^(r) as given in 
Eq. (O, Rectj^^ is the rectangular pulse function, 
which yields unity when its argument is within the first 
Brillouin zone defined by the reciprocal-lattice vectors, 
{fei, &2, ^3}, and zero otherwise. 

As will be supported by our following results, concomi- 
tant with the Stark redshift of the single-particle energies, 
the segregation of the electron and hole wave functions 
gives rise to a blueshift that partially negates the QCSE. 
To account for this effect, we include perturbatively^iS 
the so-called diagonal direct Coulomb term, Jw.cc be- 
tween the valence-state wave function, ij^yif), and the 
conduction-state wave function, ipci'r), using the expres- 



sion 

J.... = - / dV,d3,^f^M£(4p^, (7) 
J £{ri,r2) \ri - r2\ 

where, e is the electronic charge, l/e(fi,T^) is the inho- 
mogeneous inverse dielectric function for which we use 
the mask function approach of Ref. H^. This is the only, 
but by far the dominant Coulomb term that is included 
in this study. A more elaborate approach for the elec- 
trostatic as well as the nondiagonal Coulomb terms can 
be found in Ref. |4l| . Furthermore, we ignore the spin- 
orbit coupling and hence, the NC states in this work 
are doubly spin degenerate. This coupling is particu- 
larly weak in silicon with its small atomic number and as 
a matter of fact, it forms the basis for spin-based silicon 
quantum computing proposals4^ Accordingly, no spin- 
flip process is considered in the carrier relaxation follow- 
ing the optical excitation so that only spin-triplet exci- 
tons are formed for which there is no exchange Coulomb 
contribution. Even if spin flips were to be allowed, for the 
NC size ranges of this study, their contribution which de- 
cay with the third power of diameter— would become to- 
tally negligible compared to the direct Coulomb and the 
single-particle Stark energies. Obviously, future studies 
can avoid some of these simplifications in this work. 



III. RESULTS 
A. Stark effect 

In the recent experimental demonstration of QCSE in 
Si NCsfi^ the wavelength for the peak-emitted intensity 
occurs at 780 nm. Based on our prior theoretical study 
the corresponding diameter of the NC that matches with 
this optical gap is extracted as 5.6 nm. For this size of 
a NC, we display in Fig. [2] the evolution of the single- 
particle states with applied Stark field. This clearly re- 
veals that valence states are more prone to Stark shifts 
which was also observed to be the case in InP quan- 
tum dots.-iS Indeed, Fig. [3] vividly demonstrates that 
the highest-occupied molecular orbital (HOMO) wave- 
function distribution is significantly shifted by the Stark 
field and gets spatially squeezed between the high Stark 
field and the spherical NC interface. On the other hand, 
the lowest-unoccupied molecular orbital (LUMO) state 
encounters only a slight displacement, which is in the 
opposite direction with respect to HOMO as expected. 
According to stronger confinement of the valence states 
under the electric field, the interlevel separations become 
wider than those in the conduction states, as can be 
checked from Fig. [51 However, it needs to be reminded 
that the size quantization energy of the electrons is much 
larger than holes in Si NCs4i 

In Fig. m we compare the experimental Stark redshift 
datai^ at 30 K with our present theoretical results. To 
correlate with this PL experiment and account for the 
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FIG. 2. (Color online) Single-particle energy levels of a 5.6 nm 
diameter Si NC for different internal NC electric fields. 



thermal effects as well as the brightness of each excited- 
state recombination, we use the following Boltzmann 
factor-averaged and oscillator strength-averaged radia- 
tive recombination (i.e., emission) energy 



J cv 



emission 



T e 



J c 



(8) 



here, E^v = Ec — E^, E^h are the conduction (c) to 
valence (v) state and LUMO to HOMO transition en- 
ergies, respectively; f^v is the Cartesian-averaged oscil- 
lator strength of the transition^ and /3 = l/{kT) with 
k being the Boltzmann constant. We apply exactly the 
same averaging on the direct diagonal Coulomb energy 
by replacing the first Ecv term in the numerator with 
Jvv,cc- The necessity for the Coulomb term is justified 
by the excellent agreement with the experimental data 
in Fig. 131 whereas without it (i.e., at the single-particle 
level) Stark shifts become significantly overestimated. 
Since our model does not incorporate any size inhomo- 
geneity, interface states or the strain effects, its success 
also supports the atomistic quantum confinement frame- 




FIG. 3. (Color online) HOMO (upper row) and LUMO (lower 
row) wave function isosurface profiles of a 5.6 nm diameter 
Si NC under no (left column) and 0.6 MV/cm (right column) 
internal electric fields. The opposite signs of the wave function 
are represented by blue (dark) and red (light) colors. The 
electric field is horizontally directed from right to left. 



work as the main source of the luminescence in these 
particular Si NC samples. Furthermore, we note that as 
in the experimental work,^* we do not observe a dipolar 
term that gives rise to a linear dependence to the electric 
field. 

The inset of Fig. 2] shows the evolution of the va- 
lence and conduction single-particle states of a 5.6 nm- 
diameter Si NC with respect to dc electric field. It in- 
dicates that there exists a level crossing between the 
HOMO and the HOMO-1 states around an internal NC 
field of 150 kV/cm. 

Figure [5] contrasts the experimental PL peak intensity 
with the theoretical radiative recombination rate, both 
as a function of applied electric field. The nonmonotonic 
behavior of the experimental data as well as the peak in- 
tensity appearing around 0.5 MV/cm and the subsequent 
fall off beyond are all reproduced by the theory. However, 
for small Stark fields, the 300 K emission comes out to 
be stronger in the theoretical estimation. This may be 
due to thermally activated nonradiative processes that 
degrade the emission rate at higher temperatures. This 
deviation shows that further work is required to properly 
account for all thermal aspects of this problem. 



B. Polarizability of Si NCs 

The theoretical results up to now were restricted to a 
single diameter of 5.6 nm as extracted from the PL peak 
of the experimental data.— We have ignored the size dis- 
tribution of the NCs in the actual samples.— Next, we 
extract the size dependence of the polarizability of Si 
NCs, defined as AE = -{l/2)aF^ij, where AE is the 
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FIG. 4. (Color online) The comparison of theoretical and 
experimental^ Stark redshifts of a 5.6 nm diameter Si NC at 
30 K. The dotted line shows the theoretical curve without the 
direct Coulomb term included. Lines are solely for guiding 
the eye. Inset shows the single-particle Stark shifts of the 
band-edge states for the conduction (upper panel) and valence 
states (lower panel). 
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FIG. 5. (Color online) The experimental PL peak intensit)*^^ 
(upper plot) and the theoretical emission rate (lower plot) for 
Si NCs at 30 and 300 K. Lines are solely for guiding the eye. 



overall Stark shift in the energy and F^c is the elec- 
tric field inside the NC. Here, the polarizability, a, is 
taken as scalar which is in general a rank-2 tensor, how- 
ever, we have observed that the variation in Stark shift 
with respect to relative orientation of the electric field 
and the crystallographic planes of the NC or the c-axis 
of the C^y point group of the NCs^ give rise to less 
than 10 meV changes for the highest applied fields. In 
Fig. iniwe show the excitonic polarizability (i.e., with di- 
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FIG. 6. (Color online) Theoretically computed polarizability 
for Si NCs based on the single-particle and excitonic (i.e., with 
direct Coulomb term included) Stark shifts at 30 K. Solid lines 
are solely for guiding the eye. The dashed lines are cubic fits 
to the data (see text). 



rect Coulomb term included) at 30 K. For comparison 
purposes, the single-particle polarizability is also pro- 
vided where this estimate becomes highly exaggerated 
for larger diameters. Both of these curves display a non- 
monotonic behavior with respect to size which exists in 
other physical properties as well; such variations occur as 
the states, such as HOMO and LUMO, acquire different 
representations of the C31, point group for different NC 
diametersj^ As in the basic dipole polarizability, their 
overall trend can be easily fitted by a cubic dependence 
a = cDf^Q, where I?nc is the NC diameter, and with 
c = 2.436 X 10"" and 4.611 x lO'^^ C/(Vm), for the ex- 
citonic and single-particle cases, respectively; in another 
unit system they are expressed as c = 1.521 x 10"'^^ and 
2.878 X 10^5 meV/[(kV)2 cm], respectively. 



C. Intraband electroabsorption of Si NCs 

In Fig. |7| the intravalence and intraconduction state 
electroabsorption of 4 nm and 6 nm diameter Si NCs 
are shown, under and 0.6 MV/cm internal NC elec- 
tric fields. Using an anisotropic effective mass model 
and focusing only on the conduction band, de Sousa et 
al.-'^ have also modeled the intraband electroabsorption 
in Si NCs and as in this work, they obtained a blueshift. 
However, as expected from the rigidity of the conduc- 
tion states under the electric field (cf.. Figs. |2|and|3]), 
and as observed in Fig. [71 the electroabsorption effect 
can be best utilized in p-doped Si NCs. Unlike the inter- 
band transitions, we obtain a blueshift in the spectra with 
the applied electric field; observe from Fig. [2] that as the 
electric field increases, both the valence and conduction 
states individually fans out, whereas the optical band 
gap gets redshifted. For the 6 nm case, the first peaks in 
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FIG. 7. (Color online) The intraconduction and intravalence 
state electroabsorption curves for 4 nm and 6 nm diameter Si 
NC at 300 K. For amplitude comparison, all four curves are 
drawn to scale among themselves. 



the intravalence electroabsorption spectra shift close to 



38 meV under a NC field of 0.6 MV/cm; this shift reduces 
to about 14 meV for the 4 nm case. Given that bulk sili- 
con has very poor Franz-Keldysh and Kerr effect efficien- 
cies for electro-optic modulation,^ these results can be 
encouraging for the consideration of nanocrystalline Si- 
based infrared electroabsorption modulators. However, 
we should also remark that the doping of Si NCs has its 
own challenges compared to bulk."*® Finally we should re- 
mark that as mentioned within the context of the Stark 
field, the surface polarization effect, also known as local 
field effect, due to dielectric discontinuity, plays a role in 
the optical absorption as well42r— However, this effect 
becomes quite insignificant for Si NCs with diameters 
larger than about 4 nm and furthermore it only affects 
the amplitude of the absorption without modifying its 
spectral profile4^ 

IV. CONCLUSION 

In conclusion, we show that the recent QCSE datai^ 
for the embedded Si NCs under a strong Stark field can 
essentially be explained very well with an atomistic pseu- 
dopotential model which otherwise does not incorporate 
any size inhomogeneity, interface states, or the strain 
effects. In this context, the importance of the direct 
Coulomb interaction is demonstrated, and a simple ex- 
pression for the Si NC polarizability is extracted. Finally, 
in compliance with the fact that the valence states display 
much more Stark shift, it is shown that intravalence band 
electroabsorption enjoys wider voltage tunability which 
can be harnessed for Si-based electroabsorption modula- 
tor intentions. 
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